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Method and apparatus for filtering out an E«Bgfall FSlflO IS SAY 2006 
Description 

5 The present invention is directed to a method according to the preamble to claim 1 , 
device according to the preamble to claim 7 and a computer program product ac- 
cording to the preamble to claim 13. 

Sensing of EMG-signals in a patient's diaphragm by placing a catheter with a num- 
10 ber of electrodes in esophagus is a known technique, which is described in, among 
others, US 5,671,752. The EMG-signals received can be used in connection with 
mechanical ventilation of patients, which among others is described in US 5, 
820,560 and WO 98/48877. 

15 Sensing of EMG-signals from the diaphragm can even take place outside the body 
with electrodes placed on the patient, such as described in e.g. US 4,248,240. 

If a catheter with electrodes is guided down in the esophagus the electrodes lie on 
both sides of the diaphragm and at different distances therefrom. Each electrode's 
20 position relative to the diaphragm is normally not known and furthermore can vary 
when the patient breathes or moves in another way. 

Since the EMG-signal from the diaphragm is relatively weak, in particular compared 
to interferences from EKG, a continuous desire is to in the best way attain the high- 
25 est quality possible signal handling of the raw signal which the sensors detect. This 
is evident even in WO 01/03579. In WO 01/03579 it is assumed that the electrodes' 
location in relation to the center of the diaphragm is known. Then the electrodes are 
weighted based on location and symmetry, in which the EKG signal is taken into ac- 
count in a traditional way. 
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Previously known methods for compensating for the EKG signal include, among 
others, using a band pass filter which filters out the fi*equencies where the EKG sig- 
nal normally appears. It is also known to measure the EKG signal separately and 
then remove an equivalent signal fi-om the measured EMG signal. 

5 

None of these methods consider the actual disturbance the EKG signal creates in a 
particular measuring situation in a particular patient. This disturbance also varies 
with time. 

10 The known methods also fail to consider that the electrodes location relative to the 
diaphragm often varies during a single measurement. 

It is a goal of the present invention to provide a method and a device for measuring 
EMG signals which takes into consideration that the electrodes location varies dur- 
1 5 ing a measurement. 

This goal is achieved in accordance with the invention by the method above includ- 
ing the characteristics in the characterizing part of the independent claim 1 and the 
device above including the characteristics in the characterizing part of the independ- 

y 

20 ent claim 6. The goal is also achieved with a computer program product comprising 
the elements in the characterized part of the independent claim 10. 

(, 

Advantageous designs and further developments of the method are evident from the 
independent claims. 

25 

By selectively estimating the signal-to-noise ratio for each signal channel and ar- 
ranging significance factors depending upon the signal-to-noise ratio, the channels 
with the highest quality can be used more effectively for determining the EMG- 
signal. 
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Since EKG-signals constitute a main disturbance source, the determination can in 
large part be done in relation to these signals. The signal-to-noise ratio can then be 
determined as the quotient between an estimated EKG-activity and an estimated 
EMG-activity. The quotient can be determined so that a high estimated EMG- 
5 activity is always awarded, e.g. by applying the equation 

where R represents the EMG-activity and S represents the EKG-activity. 

An example of the method and the device according to the invention is described in 
10 connection with the figures. In this connection, there is shown in: 

Fig. 1 a device according to the invention connected to a patient who is undergoing 
respiratory therapy. 

Fig. 2 a first example of electrode coupling for receiving a raw signal. 
Fig. 3 a second example of electrode couplings for receiving a raw signal, 
15 Fig. 4 a flowchart describing the steps which are included in the signal treatment 
with the method according to the invention. 

Fig. 5 a flowchart showing how an ECG-activity is estimated with the method ac- 
cording to the invention. 

Fig. 6 a flowchart showing how an EMG-activity is estimated with the method ac- 
20 cording to the invention, and 

Fig. 7 and Fig. 8 flowcharts which show how respectively an ECG activity and an 
EMG activity are estimated according to the method according to an alternative em- 
bodiment of the invention 

Fig. 9 a flowchart showing how a weighing of signals on different signal channels 
25 can be obtained with the method according to the invention. 

To clarify, the concept of signal and raw signal are defined to comprise n-1 signals 
from the n electrodes/sensors, and to the n-1 signals there are n-1 channels. 
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Fig. 1 shows a device 2 for filtering and analysis of EMG-signals according to the 
invention- The device 2 can in a known way be connected to a patient 4 via a cathe- 
ter 6 with a number of electrodes 8A, 8B, 8C, 8D in one end (four electrodes are 
shown, but the number can be larger or smaller). By placing the catheter 6 in 
5 esophagus (gullet) of the patient 4, the electrodes 8A, 83, 8C, 8D can be placed in 
different locations in the diaphragm 10 (the size of which has been exaggerated in 
the figure to indicate the relative placement of the electrodes 8A, 8B, 8C, 8D). In an 
analysis unit 12 in the device 2, filtering and analysis of the raw signal fi-om the 
catheter 6 is done to extract the highest quality EMG-signal possible, 

10 

In this connection, the raw signal can be received in many different ways. Figs. 2 
and 3 show two ways. From Fig. 2 it is evident that the electrodes 8A, 8B, 8C, 8D 
are coupled together in pairs via three couplers 14A, 14B, 14C and in that way give 
rise to a three-channel raw signal (with e.g. nine electrodes, eight channels are re- 
15 ceived). 

In Fig. 3 an example is shown where the respective electrodes 8A, 8B, 8C, 8D are 
connected to a reference electrode 16 (which e.g. can be grounded) via four couplers 
14D, 14E, 14F, 14G. This gives rise to a four-channel raw signal (for eight channels 
20 in this arrangement, eight electrodes and one reference are required). 

Naturally an arbitrary number of electrodes be used and for n electrodes n-1 signals ( 
are allowed and thereby n-1 channels. 

25 More information regarding the catheter, the sensors and the entire process to cap- 
ture raw signals from the diaphragm via esophagus can be found in e.g. US 
5,671,752 and WO 01/03579. As already noted, electrodes connected outside the 
body can even be used instead of completely non-mvasive receiving of EMG- 
signals. 
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The patient 4 can also be connected in conventional ways to a ventilator system 18, 
which in turn can be connected to the device 2. The respiratory therapy given via the 
ventilator system 18 in that way can be influenced by the EMG-signal, which is ex- 
tracted from the raw signal from the diaphragm 10. This influence can be done in 
5 many different ways, of which some are described in US 5,820,560 and WO 
99/43374. 

The present invention is directed to the device 2 and, to be precise, the analysis unit 
12. The analysis unit 12 filters EMG-signals out of the raw signal from the catheter 
10 6. In that connection a number of signal channels are used, as noted above. 

To receive the highest quality possible in the EMG-signals, the raw signal contains 
in addition to EMG even EKG, alternating current noise, noise, movement artefacts 
and other low frequency disturbemces, filtering in the analysis unit 12 is done ac- 
15 cording to the method described below, which can be performed analog, digital or as 
a combination of these and be realized in hardware, software or in a combination of 
these. 

The method is described in connection with the flowcharts in Figs. 4, 5, 6, 7, 8 and 
20 9. 

Fig. 4 shows a flowchart which describes the overall signal handling in the analysis 
unit 12. The general handling is as follows. 

25 The raw signals from the sensors are input by a (multi-channel) input 22 and first 
pass a first high pass filter 24. The purpose of this is to filter away movement arti- 
facts and other low-frequency interferences. The breaking frequency should be 
lower than 10 Hz. 
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In the next step, the signals are filtered in a non-linear low pass filter 26. The pur- 
pose of this is to smooth the signal when it has a high amplitude, which typically 
happens with the existence of EKG-signals. The breaking frequency should lie 
within an interval of about 50-700 Hz. Where in the interval the breaking frequency 
5 lies depends on the energy or the amplitude of the signal. Higher amplitude results 
in lower breaking fi-equency. The non-linear low pass filter 26 should even have a 
dynamic component, namely that the breaking frequency changes with a time con- 
stant. 

10 Subsequently, the signals pass a second high pass filter 28. The purpose of the sec- 
ond high pass filter 28 is to select the frequency interval where the EMG-signal lies. 
The breaking frequency is therefore chosen to the lower regions of the bandwidth of 
the EMG-signal which is about 100 Hz. 

15 The block 26 and 28 can be replaced with an adaptive band pass filter which is de- 
scribed in our parallel application 0303061-6, in connection to figure 6 in that appli- 
cation. 



In the next step, nearby channels are differentiated from each other in a differenti- ^ 
20 ator 30. The purpose of this is to remove ringing in the filter and is based on the as- 
sumption that nearby channels are correlated with reference to common-mode inter- 
ferences. I 

Subsequently, the energy content of the signal is determined in an RMS-former 32 
25 (Root Mean Square). To reduce any remaining spikes in EKG-signals (foremost re- 
lated to the QRS-complex in the EKG-signal), the derivative of the signal fi-om the 
RMS-former 32 is limited in a rate limit block 34. 
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the respective channels by a weight factor (see below), summing and normalizing 
the signals. In this connection the weight factor can be squared to more selectively 
promote the channels with good signal-to-noise ratio. In principle, the sununator can 
be seen as a channel selector where the channels which have the highest weight fac- 
5 tor are chosen out while the channels with poorer SNR can be allocated the value 0 
in extreme cases. 

The weight factor for the respective channels is determined as follows. 

10 After the first high pass filter 24, the signals in the channels are also passed to an 
EKG-detector 38. The purpose of the EKG-detector 38 is to establish the presence 
of the EKG-signals in the channels. If a single channel indicates the presence of an 
EKG-signal, an estimated EKG-activity is calculated for all of the channels. The de- 
termination of the estimated EKG-activity is done in a first calculation block 40. 

15 The output signal from the EKG-detector 38 is even led to a second calculation 
block 42. More on these below. 

In a corresponding way as described above, the presence of EMG-signals can also 
be established. This is done in an EMG-detector 44, which is fed with the signal 
20 from the second high pass filter 28. 

The output signals from the EMG-detector 44 are led to the second calculation 
block 42. 

25 In addition to the signals mentioned above, the first calculation block 40 and the 
second calculation block 42 have a further input signal, namely the signal after the 
rate limit block 34 
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The EKG detector 38 and the EMG detector 44 can be designed in different ways. 
In one formulation the EKG detector 38 is designed to, for each channel, that is, for 
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each signal which is received from an electrode to detect if the EKG signal exceeds 
a limit value which is defined for the EKG signals. To make this comparison the 
raw signal is filtered in a band pass filter to take away the relevant frequency band 
for an EKG signal and the output signal from the band pass filter is compared with 
5 the set threshold. If the output signal is higher than the threshold the EKG signal is 
considered present. 

In an alternative embodiment both the detectors 38, 44 are designed so that a first 
probability function Pekg is determined for each channel. Pekg reflects the probabil- 
10 ity that an EKG signal is present in the signal. A second probability function Pemg 
1- Pekg reflects the probability that an EMG signal is present in the signal. 

The probability functions are calculated from a frequency analysis for the respective 
EKG and EMG signals. 

15 

The function of the first calculation block 40 is evident from Fig. 5. The signal from 
the EKG-detector 38 goes into a first decision block 46. Here it is established 
whether an EKG-signal is present in one of the channels (output yes) or not (output 
no). If there is no EKG-activity, the activity level is set to 0 in block 48 (estimated 

20 EKG-activity = 0). If an EKG-signal is present in one of the channels, the estimated 
EKG-activity S is calculated for all the channels, which is done via a low pass filter 
50, which also receives the filtered signal from the rate limit block 34. In the low 
pass filter 50, the D.C. tension level for each channel is in principle determined, the 
filtered signal from rate limit block 34, (a breaking frequency of a few Hz can be 

25 accepted), which in that connection represents the estimated EKG-activity S. 

The function of the second calculation block 42 is evident from Fig. 6. The signal 
from the EKG-detector 38 goes into a second decision block 52. If an EKG-signal 
exists (output yes), the estimated EMG-activity is set to 0 in block 54. If no EKG- 
30 signal exists (output no), it is investigated whether there exists some EMG-signal 
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(from the EMG-detector 44) in a third decision block 56. If no EMG-signal exists 
(output no), the activity is set to 0 in the block 54. If there exist EMG-signals, the 
estimated EMG-activity R is calculated by passing the signal, the filtered signal 
from Rate limit block 34, through a low pass filter 58. 

5 

Another embodiment of the blocks 40 and 42 are shown in figures 7 and 8, respec- 
tively. 

Figure 7 shows an alternative embodiment of the first calculation block 40 for esti- 
10 mating the EKG signal. From the rate limit block 34 a low pass filter 50 in the first 
calculation block 40 receives the filtered raw signal. From the EKG detector 38 it 
receives as previously the estimation of the presence of EKG signal. In contrast to 
the above-described embodiment the value which is received from the EKG detector 
is a probability function Pecg with a value between 0 and 1. The output signal from 
15 the calculation block 40 is an estimation of the effect contents in the EKG signal, 
which depends on the probability function. 

Figure 8 shows an alternative embodiment of the second calculation block 42 for 
estimating the EMG signal. The output signal from the EKG detector 38, Pecg is 

20 even fed to a second calculation block 42 where the probability 1-Pecg is fed to a 
low pass filter 58. The low pass filter 58 even receives the filtered signal from the 
rate limit block 34. Before the filtered signal is fed to the low pass filter the noise 
level in the signal is estimated in an estimating block 59 and drawn from the signal. 
The output signal from the calculation block 42 is an estimation of the effect con- 

25 tents in the EMG signal, which depends on the mentioned probability fimction. 

As shown, the embodiments shown in figures 5 and 6 can be seen as a special case 
of that shown in figures 7 and 8, where the probability Pecg can assume the value 0 
or 1. 
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The estimated EMG-activity R and the estimated EKG-activity S are transferred to a 
SN-block 60 for determining a signal-to-noise ratio for each channel between the 
estimated EMG-activity R and the estimated EKG-activity S. The signal-to-noise 
ratio can be determined in different ways. One way to determine a signal-to-noise 
ratio T is through the quotient 



where R is the estimated EMG-activity and S is the estimated EKG-activity. 

To promote channels with high estimated EMG-activity (regardless of interfer- 
ences), it is advantageous to instead use the quotient 

R + S 

Instead of R^, R" can naturally be used, where n suitably is considerably larger or 
equal to 1. 

The signal-to-noise ratio is transferred to a weight-factor block 62, where a weigh- 
ing is determined for each channel. The determinations in the weight-factor block 
62 are apparent in Fig. 7. A time block 64 counts time intervals t, for example a few 
seconds long. During the respective time interval t, the maximum signal Tl is de- 
termined for the signal-to-noise ratio T from the SN-block 60 in a maximizing block 
66. This maximum signal Tl is then filtered in a low pass filter 68. The filtered sig- 
nal T2 then represents the base for the determination of a weight factor for each 
channel in a calculation block 70. In this example, the weight factor is set to 

T2 
max 72 

where maxT2 is the maximum T2 for all the channels. In other words, the channels 
are normalized to the strongest signal-to-noise ratio of all the channels, such that the 
weight factor for the respective channel receives a value between 0 and 1 . 
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The determined weight factors are then transferred to the summator 36 (Fig. 4), 
where handling is done as previously described. Finally, the signal can be smoothed 
in a low pass filter 72 and put out as a ready EMG-signal. 

5 Preferably, a computer program is arranged in the analysis unit 12, which directs the 
function of the different blocks in the analysis unit according to the above. 

In the method according to the present invention, the determination of weight fac- 
tors for the respective signal channels constitute the main portion. Still, it can again 
10 be pointed out that the numerical values given above, e.g. breaking frequencies and 
multiplication factors, are only given as examples and in no way exclude other val- 
ues. Similarly, certain portions of the signal handling can be obtained in other ways 
without deviating from the essence of the invention. 



